Prom non-Brownian Functionals to a Fractional Schrodinger Equation 
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We derive backward and forward fractional Schrodinger type of equations for the distribution 
of functionals of the path of a particle undergoing anomalous diffusion. Fractional substantial 
derivatives introduced by Friedrich and co-workers [PRL 96, 230601 (2006)] provide the correct 
fractional framework for the problem at hand. In the limit of normal diffusion we recover the 
Feynman-Kac treatment of Brownian functionals. For applications, we calculate the distribution of 
occupation times in half space and show how statistics of anomalous functionals is related to weak 
ergodicity breaking. 
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Brownian functionals have many applications in 
physics, and hence are well investigated [l|, 0- With 
the path of a Brownian particle x(t) in the time inter- 
val (0,i), we define the functional as A = J Q U[x(t)]dt, 
where U(x) is some prescribed function. Since x(t) is 
a random path, A is a random variable. As discussed 
in a recent review , Brownian functionals model many 
phenomena, such as: fluctuating interfaces (U(x) = x 2 ), 
statistics of occupation times (U(x) is equal to the step 
function 9(x)), Obukhov's model of advection of particles 
in turbulent flow 0, 0| (U(x) = x and the velocity field 
is modeled with Brownian motion) , and finance of stock 
prices (U(x) = exp(— /3a;)) |2j, to name only a few exam- 
ples. Kac used Feynman's path integral method to obtain 
the (imaginary time) Schrodinger equation for the distri- 
bution function of A P, . The celebrated Feynman-Kac 
formula is based on the assumption that the diffusion is 
normal Brownian motion. However, we know today that 
in a vast number of applications in physics the underlying 
processes exhibit anomalous diffusion 0,0,[3|- Hence, the 
whole machinery of treating Brownian functionals should 
be modified for anomalous cases. 

Several authors have recently investigated specific 
functionals of anomalous processes. Examples include 
local and occupation times (U(x) = S(x), U(x) — 9{x), 
respectively) in random walks in the Sinai model 0], a 
fractional Obukhov model (U(x) = x) 0, and statistics 
of occupation times (U(x) — 0(x)) of a particle under- 
going sub-diffusive process in half space 0] (see more 
details below). In this manuscript we obtain a rather 
general framework for anomalous functionals, thus gener- 
alizing the Feynman-Kac formalism to an important class 
of stochastic processes. We use the continuous time ran- 
dom walk model [loj ] . which is a well established model 
for anomalous diffusion with many applications in physics 
[fil, 0, E3] • Specifically, we will soon derive forward and 
backward fractional Schrodinger equations, where the or- 
dinary time derivative is replaced by a fractional substan- 
tial derivative [HI], an d the function U(x) plays the role 
of a potential. A few applications are then worked out. 

CTRW model. We consider a random walk on a one 



dimensional lattice with lattice spacing a. Jumps are 
to nearest neighbors only and with equal probability of 
jumping left or right. Waiting times between jump events 
are independent identically distributed random variables 
with a probability density function (PDF) ip(r). Thus, 
the particle waits on lattice point xo (the starting point) 
for time r drawn from ip(r), then jumps with proba- 
bility 1/2 to xq + a or xq — a, and then the process is 
renewed. The main interest of this Letter is the case 
iP(t) ~ B Q T-( 1+Q! V|r(-a)|, when < a < 1. In this 
case, the waiting time is moment-less ((t) = oo) and the 
diffusion becomes anomalous. Values of a for a large 
number of systems and models are given in 0, 0, [13] • 

CTRW Functionals. Let G(x,A,t) be the joint proba- 
bility density function of finding the particle on (x, A) at 
time t. Denote the time the particle performed the last 
jump in the sequence as t—r. According to the model the 
particle is on (x, A) at time t if it was on [x, A — tU (x)) 
at time t — r, immediately after the last jump was made. 
Let Q n (x, A, t)dxdA be the probability per unit time, to 
arrive in [(ar, x + da:), (A, A + dA)] after n jumps. Then 



G(x,A,t) 



W (r) V Q n [x, A - tU{x), t - t] dr, 



n=0 



(1) 

where W(t) = 1 — J Q ^)(r')dr' is the probability for not 
moving in time interval (t—r, t). The summation over n is 
over the random number of jumps made. The probability 
of arriving in (x,A) after n + 1 jumps Q n+ \{x, A,t), is 
recursively related to Q n (x, A,t) according to 

Q n +i(x, A, t) = f Q %1>(t) I {Q n [x - a, A - tU (x - a) ,t - 
+Q n [x + a, A — tU (x + a) , t — r] } dr. 

. (2) 

To derive this equation, we notice that to reach point x 
after n + 1 steps one has to be after n steps, in one of the 
nearest neighbors (x — a or x + a) with probability 1/2 
for each event. Hence to reach A, one has to jump from 
A — tU (x ± a), where here r, the time interval between 
jumps, is randomly distributed with the PDF ip(r). 

In what follows we consider U(x) > 0, namely func- 
tionals with positive support, where it is natural to in- 
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troduce the Laplace transform A — > p [12| . It is easy to 
see that 

Q n [x,A-U(x)T,t-T]e- pA dA = e-P u ^ r Q n (x,p,t-T), 

(3) 

where along this work we use the convention that the 
variables in the parenthesis define the space we are work- 
ing in, and thus Q n (x,p, t) is the Laplace A — > p trans- 
form of Q n {x, A, t). We now Laplace transform Eq. ([I]) 
with respect to time, t — ► s, using the convolution theo- 
rem and Eq. §5§ to find 

&{x,p,s) = , Z¥77Za Qn{x,p,s), (4) 



n=0 



S + pU (x) 



where "ijj{s) — tp(t) exp(— st)dt is the Laplace trans- 
form of the waiting time PDF. Fourier transform, x — > k, 
of Eq. (4} yields 

G(k,p, S ) = f2 1 -^ [S ffi ( 7f )] Q»(fc, P , (5) 



since G(k,p = 0, s) is the Fourier-Laplace transform of 
the PDF of finding the particle on x at time t. 

Derivation of forward fractional Schrddinger equation. 
As mentioned in the introduction we assume that the un- 
derlying process is described by power law waiting times 
with diverging first moment. For that case we have the 
small s expansion 



ifj(s) = 1 - B a s a + 



(9) 



As well known in this case, the diffusion is anomalous 
with {x 2 {t)\ = 2K a t a /T(\ + a), and K a = a 2 /2B a (units 
mt 2 /sec a ) til- Substituting Eq. © in Eq. © we find 
the long wavelength k — > and long time s — * limit 



G(k,p, s) 



K a k 2 +[s+pU(~i^)]' 
(10) 

Rearranging the expression in the last equation we find 



where the well known Fourier transformation x — > — i-^ 
was used. To complete this part of the derivation we 
calculate Q n (k,p, s) using Eq. © and find the iteration 
rule 



i + pt/(-»|)] G(k,p,s)-1 = 
K a k 2 [a+pU(-i-fa)] 1 - a G(k,p,s). 



(11) 



Qn+i (k,p,s) = cos(ka)4> 



Q n (k,p, a). 

(6) 

The initial input for n = is Qo(k,p,s) = 1, since 
Qo(x,A,t) — 6(x)S(A)S(t), namely the particle is on 
x = and A = at time t = 0, when the process be- 
gins. The iteration rule gives Qi(k,p, s) = cos(fca)V>(s + 
pU(0)). This makes perfect sense since, just after the 
first jump we have A — U (0)r and x is either on +a 
or —a (which comes from the inverse Fourier transform 
of cos(fea)). Notice that the order of the operators in 
Eq. © is important, since cos(ka) does not commute 
with $ [s +pU(— i-§^)\ ■ This order of operators is natu- 
ral, since in CTRW we first wait and then make a jump. 
Summing Eq. ([5]) over the number of jumps n, and using 
Eq. ([6]), we find the formal solution 

G(k,p,s) = 



l-^[s+pU{-i-jL)\ l 

s + P U Hm) 1 " cos(fca)^ [s + PU (-*&)] ' 

(7) 

When p = 0, we find the well known Montroll- Weiss 
equation 0, HH 



G(k,p=0,s) 



l-V(s) 



1 



1 — cos(ka)ip(s) 



(8) 



Inverting to the space- time domain s — > t and k 
find the fractional Schrodinger equation 



dG(x,p,t) d 2 ^x- an( « \ 



x we 



dt 



pU(x)G(x,p,t). 

(12) 

The operator qD\ is a fractional Riemann-Liouville 
substantial derivative, introduced by Friedrich and 
coworkers [ll|. In Laplace space, (yD\~ a — + (s + 
pU{x)) 1 ~ a . When a = 1, Eq. JTU) is the Schrodinger 
type of equation for Brownian functionals, originally de- 
rived from Feynman-Kac formula Also note that if 
p = 0, Eq. (fT2l) reduces to the expected fractional diffu- 
sion equation [7( • Using Eq. (ITS]) , it is easy to show that 
^ = (U(x)), where (U(x)) = J U(x)P(x, t)dx, and 
P(x,t) is the PDF of finding the particle on x, i.e. the 
solution of the fractional diffusion equation. Of course 
this is the expected behavior since dA/dt = U\x{tj\. 

Backward Schrddinger equation. We now derive a 
backward equation which turns out to be very useful [l5| . 
Let G Xo (A, t) be the PDF of the functional A, when the 
process starts on x . According to the CTRW model, the 
particle, after its first jump at time r, goes through ei- 
ther xq + a or xq — a. Alternatively, the particle does not 
move at all during the measurement time (0,i). Trans- 
lating this observation to an equation, we have 
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G X0 {A,t)= / dTi>(T)-{G Xo+a [A-TU(xo),t-T}+G Xo - a [A-TU(x ),t-T}} + W(t)5[A-U(x )t}, (13) 



r 



where tU{xq) is the contribution to A from the pausing 
time on x$ in the time interval (0,t). The last term 
on the right hand side of Eq. (fl~3f describes motionless 
particles, hence the Dirac delta function. Using Laplace 
transform technique similar to that used in the derivation 
of the forward equation, we find in the continuum limit, 
the backward fractional Schrodinger equation 



9G X0 (p, t) _ x _ a d 2 



2 G Xo (p, t)-pU (x Q )G Xo (j>, t). 

(14) 

Here, the fractional substantial derivative qT>\ 01 is in 
Laplace t — > s space (s + pU(xo)) 1 ^ . Notice that this 
operator appears to the left of the Laplacian d 2 /d(xo) 2 in 
Eq. |[I4"]). in contrast to the forward equation (JT^J. Since 
in Eq. (|14[) the operators depend on xo and not on x, Eq. 
(|14p is called a backward equation. Eqs. (| 1 21 14|) are the 
main equations of this manuscript since they provide a 
general framework for treating functionals of anomalous 
processes. When a = 1 both equations reduce to the 
usual Schrodinger equations found in the Feynman-Kac 
treatment of Brownian functionals 0] • 

First Illustration: Lamperti's law for distribution of 
occupation times. We consider the occupation time in 
half space: A — J Q 8[x(r)]dT, usually denoted with t + . 
This distribution was first computed by Lamperti using 
probabilistic methods [l6| (see also [§]). Substituting 
U(xo) = 0{xq) in Eq. (fl4| . we solve the backward equa- 
tion separately for xq > and xq < 0, demanding the 
continuity of the solution and its derivative. In addition 
we have hin^^oo G Xa (p, s) = (s + p)^ 1 since if Xq — > oo, 
the particle is always in x > and thus A = t. Similarly 
lim a;0 ^_ oo G Xo (p, s) = s _1 since then A = 0. At least 
in Laplace space, the solution of the fractional equation 
is easily obtained. For initial condition xq = 0, we find, 
after inversion of the solution to the time domain, the 
PDF of the scaled functional p + = t + jt i.e. the fraction 
of time spent in half of the space 



iV 



(p+)° /2 " 1 (i-p + ) 



+ W2-1 



(15) 



(l-p+) Q + (p+) Q +2(l-p+)°/ 2 (p+)°/ 2 cos ' 



Naively, one expects that the particle spends half of the 
time in x > 0. In contrast, the Lamperti PDF Eq. (fT5|) 
has two peaks at p + = 1 and p + = 0, while its minimum 
on p + = 1/2 coincides with its average (p + ) = 1/2. In 



the limit a — > we get two delta functions on p + = 1 
and p + = 0, indicating that the particle is localized in 
either x > or x < for the whole observation time. For 
a — > 1 we recover the well known arcsine law of P. Levy 

Second Illustration: Weak Ergodicity breaking. The 
technique used so far is valid for free sub-diffusion. Our 
approach can be easily extended to study functionals 
in an external binding field, to give new insights on 
the problem of ergodicity. We consider CTRW in the 
harmonic potential V(x) = muj 2 x 2 /2 as modeled by 
the fractional Fokker-Planck equation 0, [13, 12 1- (see 
[lcl ] for normal diffusion). Algorithms that generate the 
stochastic continuous trajectories were recently proposed 
[l^ . liol ]. Consider the time average x — f Q x(r)dT/t, 

hence A — J* x{t)At and U{x) = x 12|. If the process 



is ergodic, then, due to the symmetry of the harmonic 
field, the time average x — J Q x(r)dr/t, in the long time 
limit, is statistically equal to zero. For anomalous sub- 
diffusion (a < 1) it is well known that ergodicity is broken 
}2lLl22j|. With our tool box for anomalous functionals, we 
now treat the problem of fluctuations of the time average 
x. 

We use the forward Eq. {HJ for A = J * x{t)6t 
with the modification that the Laplacian of free diffu- 
sion, is replaced by the Fokker-Planck operator L[ p = 
K a (-^muj 2 x/kbT + j^tz), which will be further justified 
in a longer publication and T is the temperature. Thus, 
the modified forward equation of motion is 



dG(x,A,t) 







dt 



= L fp oT>l~ a G(x, A, t) - x—.G(x, A, t) 

(16) 

Here in Laplace t — > s space we have oT>] 01 — > (s + 
x-^) x ~ a . The same equation, after renaming of vari- 
ables, describes a weakly damped kinetic model of super- 
diffusive Levy walks [ll|, [23| . We investigate the fluctua- 
tions of the time average, namely, (x 2 ) = (A 2 )/t 2 . Using 
Eq. p6p , we derive equations of motion for the low order 
moments (x), (A), (xA), (A 2 ) and (x 2 ). These equa- 
tions are closed in the sense that they are not coupled to 
any higher order moments, due to the harmonic poten- 
tial under investigation and the choice of the functional. 
Using the relaxation time t° = kbT/(K a mu> 2 ), we find 
in Laplace space 
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2v _ 2 [(l-a) + ( S rn 



2(x 2 ) th + (sr) a (xoY 



1 + (ST)" 



(sr) a 



(17) 



where (x 2 )th = kbT/mui 2 and xq is the initial condition. Inverting to the time domain, using (a; 2 ) = (A 2 )/t 2 , we find 

{x 2 ) = (1 - a) (x 2 ) th + 2a (2(x 2 ) th - (x a ) 2 ) E a , 3 [- (t/r) a ] + 2 (1 + a) ((x ) 2 - {x 2 ) th ) E a , 3 [-2 {t/r) a ] , (18) 

where E a ^{z) = YlnLa ^"/^(an + 3) is the Mittag-Lefflcr function 24]. To derive Eq. (fH3")l we used 
/ °° exp(-s<)Vs Qi3 [-ci a ]dt = s"~ 3 /(s Q + c) ^24]. For long times we have 



(x 2 ) ~ (1 - a) (x 2 ) th + [(3a - 1) (x 2 ) th + (1 - a) (x ) 2 ] T(3 - a)' 1 Q° + ((j)^) 



(19) 



r 




FIG. 1: Fluctuations of the time average of x{t), (x 2 ) ver- 
sus t/r for a particle in a binding Harmonic field, do not 
vanish when t — > oo thus ergodicity is broken. The excep- 
tion is the case of normal diffusion a = 1, which exhibits 
ergodic behavior. The dashed straight lines are the limit 
limt^oo^ 2 } = (1 - a)(x 2 } th - We use (x 2 ) th = 1/2, x = 
and hence {x 2 (t = 0)) = (zo) 2 = and according to Eq. (fl9|) 
for a < 1/3 (a > 1/3) solution approaches the asymptotic 
limit from below (above). 



- 2 )tft(l — ce) in agreement with 
= 1 we have ergodic behavior 
(|18ll9p give, for the first time, 



hence lim t ^ 00 (a; 2 } = { 
[22l |. Only when a 
lim^oo^ 2 ) = 0. Eqs. 
the time dependence of (x 2 ) which exhibits a t~ a power 
law convergence to the long time limit. The convergence 
to the asymptotic result is extremely slow when a — > 
as might be expected. Notice that the initial condition 
xq enters in the correction term to the long time limit, 
and in this sense initial conditions decay slowly. For 
short times we have {x 2 ) ~ (xq) 2 + [4((x 2 )th — (^o) 2 ) ~ 
2a(xo) 2 ](t / r) a /r(3 + a). We demonstrate these behav- 
iors in Fig. [1] where (x 2 ) versus t is plotted for various 
values of a. 



To conclude we have found both forward and backward 
fractional Schrodinger equations describing distributions 
of functionals of widely observed sub-diffusive processes. 
When a binding force is included, the analysis of these 
functionals leads to a new kinetic approach to weak er- 
godicity breaking, in contrast to the equilibrium frame- 
work provided in [2l], 22\. While previous work consid- 
ered specific functionals, our work provides a very general 
toolbox for this new field. 
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